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Abstract. Based on the recent discoveries of free floating 
brown dwarfs we derive estimates of the local mass density 
of this population of objects. Mass density estimates from 
various surveys span the range 0.03 to 0.005 Mq pc -3 . 
These estimates are compared with the local mass den- 
sities of the other constituents of the galactic disk and, 
in particular, with the dynamically determined total local 
mass density. We argue that brown dwarfs might indeed 
contribute significantly to the local mass budget, but that 
a local mass density as high as 0.03 Mq pc~ 3 as suggested 
by Ruiz et al. (1997) is rather unlikely 
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1. Introduction 

The recent discoveries of free floating brown dwarfs by 
the Calan-ESO proper motion survey (Ruiz et al. 1997), 
the DENIS mini survey (Delfosse et al. 1997, Tinney et 
al. 1997), the BRI survey (Irwin et al. 1991, Tinney 1998) 
and the UK Schmidt deep photographic survey (Hawkins 
et al. 1998) have now firmly established the long suspected 
existence of a population of single brown dwarfs in the 
disk of the Galaxy. Among the many interesting questions 
raised by these important discoveries is the local mass 
density of such objects. Although based at present on very 
low number statistics Ruiz et al. (1997) have pointed out 
that brown dwarfs may contribute significantly to the local 
mass budget. It is the aim of this note to assess the various 
density estimates and relate them to the total local mass 
density. 



2. Estimates of the local mass density of brown 
dwarfs 

2.1. Calan-ESO proper motion survey 

The Calan-ESO proper motion survey has at present led 
to the discovery of one brown dwarf, Kelu-1 (Ruiz et 
al. 1997). The authors estimate that their effective search 
volume had a size of 23 pc 3 , which implies a number den- 
sity of 0.04 pc -3 or a mass density of 0.0028 M© pc~ 3 , if 
a typical mass of a brown dwarf of 0.065 M© is assumed. 
However, these numbers may underestimate the actual 
density. First, the survey detects only stars with proper 
motions exceeding 0.25 arcsec/yr, so that some stars with 
low space velocities might not be found. Furthermore, the 
plate limits imply that even among the very nearby brown 
dwarfs such as Kelu-1 only the bright brown dwarfs will 
be detected. Ruiz et al. (1997) and Tinney (1998) argue 
that these objects have ages less than one Gyr. In order to 
correct for both selection effects we assume that the distri- 
bution function of brown dwarfs in phase space is the same 
as that of the luminous stars. The velocity distribution is 
then described by Schwarzschild distributions, i.e. three- 
dimensional Gaussians with anisotropic axial ratios. As is 
well known from nearby stars (Wielen 1977), the veloc- 
ity dispersions of the Schwarzschild distributions increase 
with the ages of the stars. We assume the same behaviour 
for the brown dwarfs, because the dynamical evolution of 
a population of brown dwarfs is expected to be the same 
as for luminous stars. The velocity distribution of young 
brown dwarfs with ages less than 1 Gyr is then modelled 
as 
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In Eq. (1) U,V,W denote the space velocity components 
of the brown dwarfs in the direction of the galactic center, 
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the direction of galactic rotation and towards the galactic 
north pole, respectively, ctj, cry, aw are the corresponding 
velocity dispersions. They are assumed to increase with 
age as 



7 u,v,w( T ) — °U0,V0,W0 + Ctj,V,WT , 



(2) 



like the velocity dispersions of the luminous stars (Wie- 
len 1977). The axial ratios of the velocity ellipsoid have 
been chosen according to the most recent discussion of the 
kinematics of nearby stars by Jahreifi & Wielen (1997) as 
o-y : o-y : <j w = 2:1: 0.57. The diffusion coefficients 
Ctj,v,w follow the same ratios and have an absolute value 
of Cu + Cu + C w = 6 • 1CT 7 (km/s)7yr. Initial veloc- 
ity dispersions of cr^g + cr vo + a wa = 100 (km/s) 2 have 
been assumed (Wielen 1977, Jahreifi & Wielen 1997). For 
the solar motion with respect to the local standard of rest 
standard values of (Uq, Vq, Wq) — (9, 12, 7) km/s have 
been adopted. The velocity distribution / is calculated 
by a superposition of Schwarzschild distributions for each 
generation of stars weighted by the star formation rate r s f r , 
defined per surface density, which we have assumed to be 
constant with respect to time. Stars with higher vertical 
velocity dispersions settle with larger vertical scale heights 
in the galactic gravitational potential than stars with low 
vertical velocities , so that there will be a thinning-out of 
the density of such stars at the galactic midplanc. This is 
modelled in Eq. (1) by an extra crw( T )~ 1 term, takes 
care of all the normalization constants. 

The expected number density of brown dwarfs in a 
cone towards the direction of Kchi-1 (^=307.7°, 6=36.8°) 
is given by 



drr z / d d vf(v), 



(3) 



where d max is the radial length of the probing volume 
investigated by Ruiz et al. (1997). The observed number 
density of brown dwarfs in the cone, v bs, can be expressed 
similarly to Eq. (3), if the proper motion limit of /i > 0.25 
arcsec/yr of the Calan-ESO poper motion survey is taken 
into account in the integration over the velocity compo- 
nents. To implement this proper motion limit one has to 
transform the cartesian velocity components U, V, and W 
in Eq. (3) to radial and tangential velocities and from 
that to proper motions. The integrand in Eq. (3) depends 
then explicitly on the radial distance r. We set the upper 
boundary of the r-integration at d max = 10 pc. This is 
based on the observation that Kelu-1 is very similar to 
DENIS-P J1228. 2-1547 (Tinncy et al. 1997), implying an 
absolute magnitude of Mr = 19.5 mag, which is about the 
plate limit of the Calan-ESO survey (Ruiz et al. 1997). 1 



1 We note that this places Kelu-1 at a distance of 9 pc, 
which is less than the 'astrometric' distance derived by Ruiz 
et al. (1997). However, we do not follow the authors in assum- 
ing that Kelu-1 is exactly at rest at the local standard of rest, 



Integrating Eq. (3) and its modification numerically we 
find a ratio of 
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which changes to 0.28 or 0.44, if d max is increased or de- 
creased by a factor of 1.6 according to an estimated uncer- 
tainty of the distance modulus of ±1 mag, respectively. Up 
to now we have considered only bright and young brown 
dwarfs. Very likely there are also fainter and older brown 
dwarfs, which will contribute to the mass budget as well, 
although they have not yet been found in the surveys. As- 
suming again a constant star formation rate and an age 
of the galactic disk of 10 Gyrs we estimate from Eq. (1) 

v{t < lOGyr) _ y/ <J W0 + C w ; lOGyr - ct W q , . 

i/(T<lGyr) " V^wo + C w • IGyr - a W0 ~ 

Combining this with the estimate (4) we conclude that the 
contribution of brown dwarfs to the local mass density is 
10 times the density of brown dwarfs deduced straight- 
forward from the Calan-ESO survey. The same correction 
factor is given by Ruiz et al. (1997), although the authors 
do not explain in detail how they arrive at that estimate. 

2.2. DENIS mini survey 

The DENIS mini survey with a field size of 230 deg 2 
has led to the discovery of three brown dwarf candidates 
(Dclfosse et al. 1977). High resolution spectroscopy by 
Tinncy et al. (1997) has clearly confirmed the brown dwarf 
nature of DENIS-P J1228.2-1547 by detecting a Li absorp- 
tion line. However, Tinney et al. (1997) conclude that the 
other two objects must be also very cool low-mass ob- 
jects. If we adopt the absolute K magnitude of DENIS- 
P J 1228. 2- 1547 and the plate limit in the K band given 
by Delfosse et al. (1997), we find that the search volume 
of the DENIS mini survey for brown dwarfs is 162 pc 3 . 
This implies a local number density of brown dwarfs of 
0.019 pc~ 3 . In order to account for the fainter and older 
brown dwarfs, which have not yet been detected by the 
DENIS survey, the correction factor given in Eq. (5) has 
to be applied. 

2.3. BRI survey 

High resolution spectroscopy of very late type stars of 
the BRI survey (Irwin et al. 1991) by Tinney (1998) has 
led to the identification of the brown dwarf LP 944-20 = 
BRI 0337-3535. Tinney (1996, 1998) has also measured the 
parallax, proper motion and radial velocity of this brown 



so that its proper motion solely reflects the solar motion. In- 
deed, LTT5054 and LTT5122, which lie in the same direction 
as Kelu-1 and have very similar proper motions, but for which 
parallaxes are known, can be shown to be typical disk stars 
with non-zero space velocities. 
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dwarf. Its age is estimated as about 500 Myrs (Tinney 
1998). The low space velocity of 7 km/s is consistent with 
such a young age (Wielen 1977). The BRI survey cov- 
ers 1000 deg 2 and its plate limit is at 19.0 mag in the R 
band. In the NLTT catalogue (Luyten 1979) the appar- 
ent R magnitude of LP 944-20 is given as R = 17.5 mag. 
This implies that the search volume of the BRI survey 
for brown dwarfs is 99 pc 3 . The local number density of 
brown dwarfs is then according to this determination 0.01 
pc~ 3 . Since LP 944-20 is younger than the brown dwarfs 
found in the other surveys (Tinney 1998), the correction 
factor for taking into account the fainter and older brown 
dwarfs is 6.8 (cf. Eq. (5)). 

2.4- UK Schmidt survey 

Spectroscopy of very red stars of the UK Schmidt survey 
of ESO/SERC field 287 (Hawkins et al. 1998) has led to 
the identification of three brown dwarf candidates at dis- 
tances between 37 and 48 pc. The plate limit is at 23.1 
mag in the R band, and the field area is 25 deg 2 , although 
crowding reduces the effective area of the survey by as 
much as a factor of 4. Using the apparent R magnitudes 
and the parallaxes given by Hawkins et al. (1998) we esti- 
mate that the survey is complete up to a distance of 50 pc. 
This implies an effective search volume of 80 pc 3 and the 
local number density of brown dwarfs is then 0.038 pc -3 . 
Hawkins et al. (1998) assume that the ages of the brown 
dwarfs, which they have found, are about one Gyr. Thus 
the correction factor for the older brown dwarfs given in 
Eq. (5) has to be applied. 

3. Discussion and Conclusions 

The four surveys have led to number density estimates of 
brown dwarfs of 0.46 pc~ 3 , 0.076 pc~ 3 , 0.069 pc~ 3 , and 
0.15 pc -3 respectively. Assuming again a typical mass of 
a brown dwarf of 0.065 M© this corresponds to mass den- 
sities of 0.03 M pc~ 3 , 0.0049 M© pc~ 3 , 0.0045 Mq pc~ 3 
and 0.01 Mq pc -3 , respectively. We note, however, that 
despite the large differences the estimates are still statisti- 
cally consistent within two standard deviations. This can 
be shown by integrating the Poisson probability function 
with respect to the mean value at a given number of N 
realizations. The range of mean values (xj, x u ) of Pois- 
son distributions from which the N realizations have been 
drawn with a probability of 95%, for instance, is then given 

by 

/ —e- x dx = 0.025 . (6) 
Jo,x u N\ 

If N=l, (xi, x u ) = (0.24, 5.6) and in the case of N=3 
(xi, x u ) = (1.1, 8.8). This shows immediately that the 
statistical uncertainties of the density estimates of brown 
dwarfs overlap at two standard deviations. 



The high estimate of the space density of Kelu-1 type 
objects, leads us to ask if they would be detectable in deep 
star count data taken with the Space Telescope (HST). We 
first consider star counts in the Hubble Deep Field (HDF) 
(Flynn et al. 1996). No very red stars, i.e. stars which ap- 
peared only in the I-band images and not in the V-band 
images, were detected in the HDF to a limiting magnitude 
of 1=26.3. The number density implied by Kelu-1 itself is 
0.1 pc -3 , while figuring in the older brown dwarfs associ- 
ated with Kelu-1 raises the estimated number density to 
0.4 pc~ 3 . Taking the absolute magnitude of Kelu-1 like 
other brown dwarfs as Mj = 16.6 and their local density 
to be 0.4 pc~ 3 and assuming an exponential scale height 
of 300 pc, we estimate that 6 Kclu-1 type objects would 
have appeared in the 4.4 square arcminutc HDF, which is 
marginally inconsistent with none being observed, while 
only 0.5 would be expected if they have an exponential 
scale height of 100 pc, which is consistent with none ob- 
served. Stronger limits can be obtained using faint HST 
star counts in the Groth Strip (Gould et al. 1997), which 
covers 114.0 square arcminutes to a limiting I-band mag- 
nitude of / = 23.8. In this field we would expect 5 Kelu-1 
type objects for a scale height of 100 pc and 15 for a 
scale height of 300 pc, whereas only 1 very red object (V- 
I>5.0) was detected. However, this assumes that all of the 
0.4 pc -3 local space density of Kelu-type brown dwarfs 
would be as bright as Kelu-1 itself, whereas most will be 
older and fainter. For example, in a two-component toy 
model, containing a young component with a space den- 
sity of 0.1 pc -3 , scale height 100 pc and Mj = 16.6 and 
an old component with a space density of 0.3 pc -3 , scale 
height 300 pc and M/ = 18.7 we would expect circa 1 star 
from each component in the Groth strip, consistent with 
1 observed red star. We conclude that faint star counts in 
HDF and the Groth strip are unable to rule out the high 
space density of Kelu-1 type objects measured by Ruiz et 
al. (1997). 

These density estimates can be compared with mea- 
surements of the densities of the other constituents of the 
galactic disk. Jahreifi & Wielen (1997) derive from a dis- 
cussion of the most recent data of nearby stars that the 
mass density of luminous stars in the solar neighbourhood 
is 0.039 Mq pc~ 3 . Thus, if the mass density estimate of 
brown dwarfs based on the Calan-ESO survey is correct, 
brown dwarfs contribute almost as much to the local mass 
budget as the luminous stars. The mass density estimates 
of brown dwarfs can be also compared with dynamically 
determined local mass densities. Such measurements in- 
clude all gravitating matter. Fuchs & Wielen (1993) and 
Flynn & Fuchs (1994) have used samples of K dwarfs and 
K giants, respectively, to measure the local slope of the 
K z force law. Both measurements have led to a total local 
mass densities of 0.1 Mq pc -3 with an estimated uncer- 
tainty of 20%. A repetition of the measurement using im- 
proved data of K giants indicates a total local mass density 
at the lower end of this range (Flynn & Fuchs, in prepa- 
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ration). Creze et al. (1997) have used recently Hipparcos 
observations of A stars and determine a value of the total 
local mass density of 0.076 ± 0.015 M© pc~ 3 . The local 
surface density of interstellar gas in the form of cold HI 
and H2 is 6 M Q pc~ 2 (Dame 1993). If corrected for the 
presence of heavier elements and folded with a vertical 
scale height of 100 pc, this implies a local mass density 
of interstellar gas of 0.037 M pc~ 3 . The mass density 
of warm HI and ionized interstellar gas is more difficult 
to assess, but is probably only 0.003 M Q pc~ 3 (Kuijkcn 
& Gilmore 1984). We conclude from this discussion that 
there is evidence for a local mass density of brown dwarfs 
of the order of 0.01 M Q pc~ 3 , which would seem to fill the 
gap between dynamical mass determinations and the mass 
density of so far identified stellar and interstellar matter, 
whereas the mass density estimate by Ruiz et al. (1997) 
(1997) seems to be on the high side. Hopefully, further 
discoveries of brown dwarfs, such as announced by the 
2MASS survey (Kirkpatrick et al. 1998), will clarify this 
issue. 
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